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Abstract 

We first formulate the problem of optimally scheduling air traffic low with sector capacity constraints 
as a mixed integer linear program. We then use semidefinite relaxation techniques to form a convex 
relaxation of that problem. Finally, we present a randomization algorithm to further improve the quality 
of the solution. Because of the specific structure of the air traffic flow problem, the relaxation has a 
single semidefinite constraint of size dn where d is the maximum delay and n the number of flights. 

1 Introduction 

In this paper, given a schedule of flights and their routes we solve the problem of finding a new schedule 
that satisfies a list of sector capacity constraints and minimizes the total delay compared to the original 
scheduled. While the optimal routing problem with weather uncertainty and capacity constraints is essen- 
tially intractable, I NEGOSJ show that robust optimal routing of a single aircraft under weather and traffic 
uncertainty can be solved efficiently as a robust Markov dynamic programming (MDP) problem. So the 
problem we solve here should be seen as a second phase: given the routes computed using MDP, our relax- 
ation produces an optimal schedule, satisfying capacity constraints while minimizing the total delay. While 
this scheduling problem is significantly simpler than the global routing problem, it still combinatorial (and 
in fact NP-Hard), hence we need to formulate a relaxation to obtain approximate solutions efficiently. We 
can interpret this problem as particular case of job shop scheduling problem where a task corresponds to an 
aircraft occupying a sector at a certain time, with the aircraft routes correspond to lists of tasks that have to 
performed in sequence and where the sector capacities correspond to server capacities. In our case, a few 
key distinctions simplify the problem formulation. First, the different tasks are independent (one aircraft's 
route is not dependent on another). Second, for each aircraft, no time gap is allowed between two tasks (we 
assume that aircraft can't be held en-route). 

Scheduling problems are notoriously hard combinatorial problems. Classic instances include the job- 
shop scheduling problem (see [CB76| for example) or the travelling salesman problem (see | RSL77| for 
example). Semidefinite relaxations and randomization techniques have and excellent track record, which 
can be traced to ILS91I . IAli95l . |GW95| and |PW951 among others. Recently, IFMOFOll applied these 
techniques to aircraft conflict avoidance for free flight, by adjusting aircraft speed and bearing. 

In this paper, we first formulate the air traffic flow scheduling problem as a mixed integer linear pro- 
gramming problem. We then apply the lifting procedure of IGW95I to formulate a semidefinite relaxation 
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of the problem. Because of the structure of our problem here, we only have dn binary variables, where d 
is the maximum delay and n the number of flights, which means that the semidefinite relaxation has one 
semidefinite constraint of size dn + 1, allowing to scale better than classic scheduling problem relaxations. 
We then detail a randomization procedure based on that of | GW95 1 that uses the matrix solution of the 
semidefinite relaxation to further improve the quality of the solution. While IFMOFQTI focus on free flight 
and conflict avoidance, our main concern is on meeting sector capacity constraints, a key limiting factor in 
the European airspace, by adjusting delay at departure. This relative simplicity allows the algorithm to scale 
very well with the number of aircraft. Finally, because of its particular structure, we show that this relax- 
ation is amenable to first-order methods for large-scale semidefinite optimization such as those detailed in 
[HROOl and INes04l which are natural algorithms for solving large-scale problems for which a low precision 
is required. 

The paper is organized as follows. The second section defines the main scheduling problem and for- 
mulates it as a mixed integer linear program. The third section derives the semidefinite relaxation of that 
problem. In the fourth section, we show how to exploit the result of this relaxation to further improve the 
solution using a randomization technique. Finally, in section five, we present some numerical results. 



2 Problem formulation 

Suppose we are given routes for n aircraft flying across an airspace composed of m sectors with capacities 
given by C G R"*. We decompose a particular day into T periods, so that a particular flight route starting at 
time s can represented by a matrix G jjmxT ^^^^^ 

{Rj^^^ = 1 if aircraft i is in sector j at time t 
Rff'^ = if not. 

We can formulate the problem of minimizing total delay while satisfying capacity constraints as: 

minimize Y.7=i Ej=o ^ijj 
subject to Er=i E -=0 ^vR^'^'^'^ <C 
T'' x -1 

Xij G {0, 1}, i = 1, . . . , n, j = 0, . . . , d, 

in the variable Xij G jj"x('^+i) where d is the maximum delay (in units of time). Here, Xij = 1 means that 
aircraft i will be delayed by j units of time (because x is a binary variable and Ej=o ~ ^' there can only 
be a single nonzero Xij for each aircraft i). The first constraint makes sure that sector capacity constraints 
are met and the objective is the sum of aircraft delays. 



3 Semidefinite Relaxation 

In this section, we apply the lifting procedure detailed in IGW95I to form a semidefinite relaxation of 
problem Q- We can rewrite Q as a non-convex quadratic program (QP): 

minimize Er=iEi=o^iiJ 

subject to E?=i EU < C 

x -1 ^' 
^ij - Xij = 0, i = l,...,n, j = 0,...,d, 
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in the variable Xij S R"><('^+1). We can rewrite this problem as: 

minimize Y17=i Ej=o 

subject to ELi E,to x^jR^'^'+^^ < C 

S^=o^ii = l' i = l,...,n (3) 

X = vec(x) vec(2;)-^ 

diag(X) = vec(x) 

in the variables Xij E and X G xhis is a non-convex quadratic program and is computa- 

tionally hard. We now detail how to obtain a convex relaxation of problem ^ using Lagrangian duality. 

3.1 Lagrangian relaxation 

Let us start from a general nonconvex quadratically constrained quadratic program (QCQP): 

minimize x'^Pqx + qTx + 



subject to x'^PiX + q[x + < 0, i = 1, . . . ,m, 
with variable x G R", and parameters G S", G R", and rj G R. We form the Lagrangian, 

(m \ / m \ m 

-Po + X] AjPj + g'o + X] Aj^i X + ro + ^ Xin. 
1=1 / \ i=i J i=i 

To find the dual function, we minimize over x, using the general formula (see example 4.5 in LBVO^ '): 

inf x-Px + + r = I ^ - i^'^;'^' ^ ^ , G 7^(P) 
^.gR [ — oo, Otherwise. 

The dual function is then: 
xgR 

^ / m \^ f \^ f \ 

V i=i / \ i=i / \ i=i J i=i 

We can form the dual of Q, using Schur complements (cf. §A.5.5): 

maximize 7 + X^i Aj^j + tq 

(Po + E™iAiP.) ('Zo + E:^iA.g.)/2 



(4) 



subject to 



^0 (5) 



Aj > 0, i = 1, . . . , m, 

in the variable A G R*". As the dual to (|4]i, this is a convex program, it is in fact a semidefinite program 
(SDP). This SDP is called the Lagrangian relaxation of the nonconvex QCQP. It can be solved efficiently 
and gives a lower bound on the optimal value of the nonconvex QCQP. We form take the dual of program 
dSli (see I.BV04. §5.9.2]): 

minimize Tr(XPo) + (IqX + tq 

subject to Tr(XPi) + qjx + < 0, i=l,. . . ,m, 

^0, 



X X 
x^ 1 



with variable x G R", X G S" and parameters Pj G S", G R", and rj G R. 
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3.2 A relaxation of the scheduling problem 

Using the results of the previous section, we can form the Lagrangian relaxation of (|5Jl as: 



minimize Ya=i J2j=o ^ijj 
subject to Er=iE' 

1, i 



X 
vec(x)'^ 



l,...,n, 



vec(x) 



1 



(7) 



^ 



diag(X) = vec(x), 

which is a semidefinite program in the variables Xi^ G ^rix(d+i) ^ £ S"('^+^) and can be solved 
efficiently. The objective of this program is a lower bound on the global solution. 



3.3 First-order methods 

An important structural property of problem is that, because diag(X) = vec(x) and l^x = n the 
matrix variable: 

X \rec{x) 
vec(x)-^ 1 

has constant trace equal to n + 1. This means that the dual of (Q is a maximum conic eigenvalue mini- 
mization problem for which efficient first-order methods such as the spectral bundle algorithm of |HROO| 
and the optimal first-order method of [Nes041. In Section |5j we present numerical experiments using the 
SBmethod code by IHROOl on increasingly large randomly chosen problems. SBmethod solves large-scale 
dense instances for which it returns a solution {x, X). It can also solve much larger problems by exploiting 
sparsity, in which case however it only returns the optimal x, which somewhat decreases the performance 
of the randomization methods detailed in the next section. 



4 Randomization 

The Lagrangian relaxation techniques developed in ^3.1l provided lower bounds on the optimal value of the 
program in Q, but did not however give any particular hint on how to compute good feasible points. The 
semidefinite relaxation in ^ produces a positive semidefinite or covariance matrix together with the lower 
bound on the objective. In this section, we exploit this additional output to compute good approximate 
solutions with, in some cases, hard bounds on their suboptimality. 



4.1 Randomization 

In section im the original nonconvex QCQP: 

minimize x'^Pqx + q^x + ro 

subject to x'^PiX + qfx + < 0, i = 1, , 



was relaxed into: 



minimize Tr{XPo) + q^x + tq 
subject to Tr{X Pi) + qfx + < 0, 
X X 
1 



i = 1 



,m, 



, . . . ,m, 



(8) 



^ 0. 
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The last (Schur complement) constraint being equivalent to X — xx^ >z 0, if we suppose x and X are the 
solution to the relaxed program in then X — xx^ is a covariance matrix. If we pick x as a Gaussian vari- 
able with x ~ J\f{x, X — xx'^), X will solve the nonconvex QCQP in (|4]i "on average" over this distribution, 
meaning: 

minimize Ei{x^ P^x + x + rg) 

subject to 'E{x^PiX + qfx + rj) < 0, i = 1, . . . ,m, 

and a "good" feasible point can then be obtained by sampling x a sufficient number of times, then simply 
keeping the best feasible point. Of course the direct sampling technique above does not guarantee that a 
feasible point will be found. In particular, if the program includes an equality constraint, then this method 
will certainly fail. However, it is sometimes possible to directly project the random samples onto the feasible 
set. As we will see below, this is the case here 



4.2 Randomized schedules 



Suppose we have solved 



minimize Y17=i Ej=o ^ijj 
subject to Er=i E,to x^JR^'^'+^^ < C 
Xyj=o "^^i ~ ^' i = 1, ■ ■ ■ , 
X vec(x) 
vec(x)'^ 1 
diag(X) = vec(x), 



y 



to get an optimal Xij € R^^^C'^+i) and X G ^"{'i+i) ^^^j^ simply sample a Gaussian variable 
Af{x, X — xx^) and compute its projection v on 



u ~ 



|0^^nx(d+l)pj^^^jjnx(d+l) . i = l,... 

j=0 



n 



We then compute the delay for each feasible random schedule v and keep the best solution. 



5 Numerical Example 

To fix ideas, let us begin by solving a very simple scheduling example. Suppose that there are only four 
sectors and two flights. 
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Each sector has capacity one. In our format, a solution where both flights leave on time will be represented 
by a; = (1,0, 1,0) and a solution where flight 1 is delayed by one unit of time will be written as x = 
(0, 1, 1, 0). Flight one starts from sector 1 and ends in sector 4, while flight 2 starts from sector 3 and also 
ends in sector 4. It takes each flight one unit of time to cross each sector and both flights are scheduled to 
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depart at time 1 and conflict in sectors 2 and 3. The maximum delay is 1 unit of time. The problem variables 
are x G R"^ and X € S^, the linear sector capacity constraints impose: 
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X <1 



and the objective vector is given by c 
following solution: 



X 




(0, 1,0, 1). We solve Q using SEDUMI by IIStu99ll to get the 



and X 



0.50 
0.24 
0.24 
0.24 



0.24 
0.50 
0.24 
0.24 



0.24 
0.24 
0.50 
0.24 




Here, the solution is not binary valued. The reason for this is symmetry, the problem too simple and sym- 
metric in flights one and two, so the solution is a mix of the two optimal solutions. This is easily fixed 
by adding a small random perturbation to the objective (implicitly breaking the tie by giving one aircraft 
priority over another). The solution is then: 



V 



and X 




Here, the fact that X 



XX 



also shows that the relaxation is tight and that the solution x is globally optimal. 



5.1 Computing times 

Here, we generate random problems in an airspace with 50 sectors and a maximum delay of 2 units of time. 
In Figure n plot of CPU time (in minutes) versus number of aircraft for these problems. 

With n = 20 for example, the computing time is about 6 minutes, for a problem with 10^*^ possible 
schedules. While general-purpose interior point solvers such as SEDUMI by fStu991 can solve reasonably 
large problem instances. Using the fact that the semidefinite program in Q has constant trace, more spe- 
cialized first order methods (see LHROOl or IdEGJLOSI and ld'A05J for example) can be used to solve much 
larger problems. 



5.2 Randomization 

To illustrate the randomization method detailed in ^ we first solve problem © on a random problem 
with 5 sectors, 3 flights and a maximum delay of 4 units of time. The solution x to Q has an objective 
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Figure 1: Log-log plot of CPU time (in minutes) versus number of aircraft for randomly generated 
problems, left: interior point methods, right: spectral bundle method. 



value of 3, which gives a lower bound on the global optimum. We then generate 100 sample points from 
u ~ M{x, X — xx^) and compute their projection v on 

|Q^^nx(d+l)pj|^^ jjnx{d+l) . =1, i = l,...,n|. 

For those sample schedules that meet the capacity constraints, we compute the total delay. Figure |2l shows 
the distribution of these objective values. We notice that the best delay found by randomization is equal 




to 3, which matches the lower bound produced by solving problem Q. This shows that 3 is the globally 
optimum delay and that the corresponding randomized solution is optimal. 
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